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1. Introduction 

At the end of the 1990s it was realised that, contrary to previously held beliefs, 
effective localisation-delocalisation transitions can take place in one- dimensional (ID) 
and quasi-lD random models, provided that disorder is endowed with specific long-range 
correlations [IJ EJ [3j E] . In particular, the study of ID models with correlated random 
potentials was greatly spurred by the seminal paper of Izrailev and Krokhin [3j H], who 
derived an analytical formula for the localisation length in one-dimensional models with 
weak correlated disorder. The discovery triggered an intense research activity in the 
field of low- dimensional models, which are now understood to possess a much richer 
behaviour than once thought (for a review, see [5]). 

The theoretical predictions about the existence of delocalised states in ID or quasi- 
lD models with correlated disorder [6j [TJ EJ EJ H] soon found experimental confirmations 
in semiconductor superlattices [7] and microwave waveguides [H [9]. Since these early 
works, ID models with correlated disorder have found applications in an increasing 
number of physical fields, including Bose-Einstein condensates [10J [Hi [12], bilayered 
media [13] , structures with corrugated surfaces [HI EES] , and DNA modelling [T6J [17] . 

In most of these works, the system under consideration has a single self-correlated 
random potential. Recently, however, some theoretical studies have begun to explore 
more complex models, characterised by the simultaneous presence of two kinds of 
disorder. For instance, in [13] the authors considered the transport properties of generic 
bilayered structures formed by infinite arrays of alternating slabs of two different kinds. 
Since slabs of both types can have random widths, each structure is characterised not 
by one but by two random sequences which, in addition to self-correlations, can exhibit 
cross-correlations. Another model exhibiting two types of disorder is the aperiodic 
Kronig-Penney model with both compositional and structural disorder [18j [191 [20]. F°r 
this model, it was possible to work out a perturbative expression for the localisation 
length of the quantum states valid for weak disorder with any kind of self- and cross- 
correlations [18]. The theoretical results for the localisation length of electronic states 
showed that the interplay of structural and compositional disorder could enhance the 
delocalisation-localisation transition created by long-range self-correlations of each type 
of disorder. It was also shown that cross-correlations could affect the transmission 
properties of a long random sample, increasing or decreasing the transmittivity in the 
opaque windows corresponding to the localised states. 

These findings posed the question of whether these effects could be experimentally 
observed in appropriately designed random microwave waveguides. In fact, the structure 
of the electronic states in a Kronig-Penney model and the transmission of the waveguide 
with point-like scatterers are closely related problems because the Schrodinger equation 
for the former model has the same form of the Helmholtz equation for the latter. Because 
of this analogy, the theoretical results originally obtained in [21] for a Kronig-Penney 
model with purely structural disorder immediately found experimental applications in 
the analysis of the transmission properties of single- mode waveguides [H [91 [22]. It 
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was therefore natural to ask whether the theoretical predictions concerning the effects 
of cross-correlations in the aperiodic Kronig-Penney model could find an experimental 
confirmation in suitable random waveguides. The main goal of this paper is to analyse 
this question, which we answer in the affirmative. 

We considered microwave transmission of the first mode in a Q1D waveguide open 
on both ends and with no intermode coupling. The waveguide contained scatterers 
constituted by brass bars situated at random positions (structural disorder) and having 
random heights (compositional disorder). We found that cross-correlations of the 
random spacings with the heights have a sizable effect on the transmission properties of 
the waveguide in agreement with theoretical predictions. This is a remarkable finding, 
because the theoretical results for the localisation length are derived for an infinite array 
with weak disorder, whereas the experimental measures were obtained using a setup of 
finite size, with a modest number of scatterers (40 bars) and relatively strong disorder. 
The experimental data therefore show that the effect of disorder correlation is robust 
and is not washed away by other features of the apparatus (such as absorption) which 
are neglected in the simplified theoretical description. 

This paper is organised as follows. In section[2]we outline the theoretical framework. 
In section [3] we describe the setup. Our experimental findings are presented in section HI 
We summarise our findings and expose our conclusions in section [5j 

2. Theoretical considerations 

We analyse the transmission of microwaves in a rectangular waveguide with internal 
scatterers, constituted by brass bars situated at random positions and having random 
heights. If one assumes that the bars can be assimilated to point-like scatterers and 
represented as delta-barriers, the propagation of a single mode in the waveguide can 
be described using a Helmholtz equation which has the same mathematical form of the 
Schrodinger equation for the aperiodic Kronig-Penney model with delta-barriers [HI |9] 



■J' 



J2 (U + u n )5(x 



n 



ip{x) = Eip(x). (1) 



For the sake of simplicity, here and in the rest of the paper we use energy units such 
that h 2 /2m = 1 and set the lattice step as unit length. 

The model (CQ), as the associated waveguide, has a twofold randomness: on the 
one hand, the positions of the delta-barriers present random shifts a n with respect to 
the lattice sites (structural disorder). On the other hand, compositional disorder is 
introduced via the fluctuations u n of the barrier strengths around the mean value U . 
The statistical properties of model ([1]) are defined by assuming that the random variables 
u n and a n have vanishing averages and known binary correlators. In the case of weak- 
disorder no further specification of the statistical features of the model is required. The 
weak-disorder case is identified by the conditions jTHl EES EQ] 

«}<£/ 2 , (A2>£«1, and (A2)f/<1, 



Experimental effects of disorder cross-correlations 



4 



with A n = a n+ \ — a n representing the relative displacement of two contiguous barriers. 

If disorder is weak, one can obtain an analytical expression for the localisation 
length li oc of the electronic states ip of the model (CE]) [I8j [191 EQ]- Following the 
Hamiltonian map approach, one can show that, within the second-order approximation, 
the electronic states have an inverse localisation length equal to 

sin 2 k 



loc o • 2 

8 sin k 



{vttWxW + U'iAttWiiK) 

sin(k) 1 & 

- 2U— ^(u n A n )cos(K) W 3 (k) 

In (J2J) the symbols k and k represent respectively the propagation wavevector within the 
wells k = y/E and the Bloch wavevector k of the quantum particle. The two quantities 
are linked by the condition 

cos k = cos k + — sin k, (3) 

which defines the band structure of the Kronig-Penney model (JT]) in the absence of 
disorder. The functions Wi(n) in (|2J), on the other hand, are the Fourier transforms of 
the normalised binary correlators of the random variables u n and A n , i.e., 

W 1 (k) = l + 2£%^cos(2 K /) 

W 2 (k) = l + 2 f;%^>cos(2K0 (4) 

1=1 X^n) 

Equation shows that the localisation length diverges (within the second-order 
approximation) in any given energy interval with vanishing power spectra (jlj). Therefore 
an effective localisation-delocalisation transition can occur in the ID Kronig-Penney 
model ([1]) provided that the disorder exhibits the specific long-range correlations which 
make the power spectra (T4j) vanish over a continuous energy range. 

In order to obtain a specific localisation-delocalisation transition with pre-assigned 
mobility edges, one must know how to construct two random sequences {u n } and {A„} 
such that the corresponding power spectra (jlj) vanish in pre-defined energy intervals. 
A solution to this problem was presented in [20]; the basic idea is to generate two 
independent white-noise sequences which, after having been cross-correlated, can be 
endowed with the appropriate self-correlations with a filtering process. More precisely, 
one starts with two sequences {A^ 1 ^} and {A^ 2 ^} of independent identically distributed 
random variables. Two cross-correlated white noises can then be obtained by using the 
transformation 

Yjp = A«cos7? + X^smri 
= X^sinrj + Xp cos rj, 

where r] determines the degree of cross-correlation of the Y variables. The values of 
rj lie in the interval [— 7r/4, 7r/4] with 77 = ±7r/4 corresponding to the extreme cases 



Experimental effects of disorder cross- correlations 



of total positive and negative cross-correlation, while for 77 = the cross-correlations 
vanish. The cross-correlated white noises Yjp and can then acquire the required 
self-correlations by means of the convolution products 
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in which the coefficients ct n and (3 n are derived by the pre-defined power spectra W\ and 
W 2 via the equations 
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After substituting the random variables defined by (JSJ) and ([6j) in 
localisation length takes the form 



(6) 

the inverse 
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As a simple example of delocalisation transition produced by self- correlations and 
modulated by cross-correlations, one can consider the case in which both structural and 
compositional disorder have a power spectrum of the form 

7T 



Wi(k) = W 2 (k) 
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with < Ki < k 2 < it/2. Note that it is enough to consider the [0, 7r/2] interval in the 
previous definition because the power spectra (jlj) are even functions of period ir. The 
binary self-correlators corresponding to the power spectra (jSJ) decrease with a power law 
(ui +n ui) (A l+n Ai) 1 



[sin (2K 2 n) — sin (2« 1 n) 



(9) 



which shows that the random variables u n and A n have long-range self-correlations. To 
simplify further the physical picture, one can consider the case in which structural and 
compositional disorder have the same strength: 

(ul) = (Al) = a\ (10) 



When conditions flS} and (flOj) are met, the inverse localisation length (|7j) for k > 
reduces to the form 
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Figure 1. (Colour online) Inverse localisation length £r* versus Bloch vector k. The 
solid black line (2) represents the results obtained in the absence of cross-correlations 
between structural and compositional disorder, while the blue dotted line (1) and the 
red dashed line (3) respectively correspond to the cases of total positive and negative 
cross-correlations. 



with I = [k%, K2] U [tv — K2, 7r — K\\. In this case the electronic states are expected to be 
localised in the intervals [k\, K2] and \n — ^2, 7r — Ki] and extended (within the second- 
order approximation) in the complementary regions of the [0, it] interval (which is the 
left half of the first Brillouin zone). 

A numerical evaluation of the inverse localisation length confirms this conclusion, as 
can be seen from figure[TJ The data represented in figure [1] were obtained for K\ = 0.207T 
and k 2 = 0.2l7r, while the mean field and the disorder strength were set equal to 
U = 1 and J (v%) = J (A^j = 0.05. The long-range self-correlations © enhance 
the localisation of the electronic states in the two narrow regions [0.207T, 0.2i7r] and 
[0.797T, 0.807r] and delocalise the other states. Cross-correlations, on the other hand, 
modulate the localisation-delocalisation transition, making it more or less pronounced 
according to the sign of the cos ft sin(2?7) term in the right-hand side of ( II ip . 

Obviously, the localisation-delocalisation transition produced by long-range self- 
correlations has repercussions on the transmission properties of a finite aperiodic Kronig- 
Penney model sandwiched between two perfect leads. When the longitudinal size of the 
random sample is larger than the localisation length of the localised states, it is natural to 
expect that a continuum of delocalised states manifests itself in the form of a transparent 
energy window, while the regions of enhanced localisation become low-transmittivity 
gaps. As for cross-correlations, they can either enhance or reduce the opacity of the 
low-transmittivity intervals. Numerical studies confirmed this expectations [20]. They 
also showed that the effects of self- and cross-correlations can be detected in random 
samples of moderate size, as can be seen in figure [21 which represents the transmission 
coefficient T of a random Kronig-Penney model of N = 40 sites for a specific realisation 
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Figure 2. (Colour online) Transmission coefficient T versus Bloch wavevector k in 
an aperiodic Kronig-Penney model of N = 40 sites. The cases of total positive and 
negative cross-correlations are compared with the case in which cross-correlations are 
absent. The represented data were obtained for the same parameters of the model 
used in the case of figured! i.e., U = 1 and \J (u^) = y/ (A^) = 0.05. 



of the disorder. 

One can see that self-correlations create transmission gaps corresponding to the 
windows of localised states. Positive cross-correlations increase the transmittivity in 
the low-energy gap (almost removing the gap itself) and further reduce it in the high- 
energy gap, whereas for negative cross-correlations the inverse effect is observed in the 
two gaps. As a final comment on this issue we note that, when one considers transmission 
in short random samples, sample-to-sample fluctuations may prevent the transmission 
coefficient from following closely the behaviour of the theoretical localisation length 
for some disorder realisations. In the next section we will discuss how the effects of 
cross-correlations appear in our random waveguide. 

3. Experimental setup and measurement technique 

We describe the experimental setup and the measurement technique in three steps. We 
first present the apparatus and explain how we measured the mode-resolved microwave 
transport through the waveguide. Then we explain how we realised the random cross- 
correlated lattice inside the waveguide. Finally we discuss how the experimental 
band structure was related to that of the theoretical model ([!]) thereby establishing 
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w = 100mm 



Figure 3. (Colour online) Schematic picture of the experimental setup with top plate 
removed (not to scale). Brass bar inlays can be freely placed inside the waveguide. 
The number of bars varies according to the spacing d„ . Antennas shown here with 
parts of the top plate can be moved separately along the z-direction. The distance 
along the x-axis between the antennas is L = 1.50 m, the total waveguide length is 
2.38 m. Absorbers at the ends of the waveguide are omitted. 

a functional relation between the wavevector k x and the Bloch vector k. 
3.1. Microwave measurement 

For our experiment we used a variant of the setup employed in [23] to analyse transport 
in quasi-lD structures. The setup consists of a rectangular waveguide. The top plate 
can be lifted so that scatterers can be placed at will inside the waveguide. Figure |3] shows 
the waveguide with the top plate removed and four bars placed inside as an example of 
scattering structure. 

Two antennas are set at positions x in = and x out = 1.5 m, with one serving 
as an emitter and the other as a receiver. The distance between the antennas, 
L = 2 ut — ^in = 1.5 m represents the effective length of the waveguide, whose width 
and height are w = 100 mm and h = 8 mm. The antennas are plugged into the 
waveguide through mobile slides so that they can be moved freely from to w along the 
transversal z-direction. Microwave absorbers are placed between each antenna and the 
corresponding end of the waveguide to avoid back-reflection from the open end. The 
total length of the apparatus is 2.38 m. 

A microwave vector network analyser allowed us to measure the transmission 
coefficient of the waveguide from the input antenna at (x in , z in ) to the output antenna 
at (x ou t,Zout), i- e -5 t ne scattering matrix S(xi n , z- m ; x ont , z out ). Because the antennas are 
fixed in the ^-direction, in what follows we write S(zi n , z out ) as a short-hand notation 
for S(0,z- m ;L,z olIb ). We measured the waveguide transmission coefficient as a function 
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of the frequency v for the lowest-frequency TE component of the electromagnetic field. 
In fact, for a given frequency v there are several propagating waves with wavenumber 
| A; | = 2nv/c. Due to the confinement along the width of the waveguide, the ^-component 
of the wavevector takes discrete values, k^> = irn/w, corresponding to propagating 
modes with longitudinal wavenumbers 

= ^k 2 - {ttu/w) 2 

with n = 1,2, ... , N w . The total number N w of propagating modes is determined by 
the integer part [[■ ■ •]] of the mode parameter kw/ir: 

N w = [[kw/n]]. 

To separate the different modes we measure S(z- m , z ont ) at several positions z- m and 
z ut distributed over the whole width w. The transmission between the initial point 
(x = 0,yi) and the final point (x = L,y-i) can be decomposed into all contributing 
modes 

S {z- m , z out ) = > S nm sm sm . (12) 

~ V w J V w J 

From ffl2|) one can obtain the scattering matrix in mode representation with a twofold 
sine transform. The matrix element S nm denotes the scattering from the m-th to the 
n-th mode. 

If all propagating modes were considered, the waveguide would correspond to a 
quasi-lD open strip [23], while the transmission coefficient for the first mode would be 
related to the scattering matrix via the identity 

n w 

T\ — ^ I Sin 1 • 

n=l 

However, we found that, as long as the scatterers fill homogenously the waveguide along 
the z direction, the individual modes propagate independently of each other through 
the waveguide. Therefore we chose bars whose length perfectly matched the waveguide 
width keeping a translation invariance in the y direction. In this way, as in [23], we 
obtained a waveguide in which the mode-mode couplings were usually negligible. This 
made the scattering matrix S mn practically diagonal (so that T\ ~ |Sn| 2 ). The absence 
of intermodal scattering allowed us to restrict our measurements to the first propagating 
mode: we were thus able to put the experimental setup in correspondence with the 
strictly ID theoretical model (JTJ. Note that the equivalence of the waveguide to a ID 
device did not require that the setup itself be strictly ID. 

3.2. Correlated random scattering structures 

We used brass bars as scatterers in the waveguide. As explained in section 13.11 the 
requirement that the setup should behave as a ID device defined the length of the 
scattering bars. The bar width b = 8 mm was chosen as in [23], because it is a good 
trade-off between small width and sufficient scattering strength. The width of the bars 



Experimental effects of disorder cross-correlations 



10 



18- 



-3)0- 

■ i — i 

CD 

a 



l.illl.illl.illi.illi.illi.illi.ilLlll 


lull 


hi 1 


Li 1 


ll.ll 


ll.ll 


ll.ll 


III 1 


■■■ill 


1. .ill. ill 


l.il 


ll.ll 


lull 


hill 


liil 


lini 







Position [mml 



800 



Positive 
cross-correlation 

No 

cross-correlation 

Negative 
cross-correlation 



Figure 4. (Colour online) Profile of the scattering arrangement within the waveguide 
for the cases of absent (top), positive (middle), and negative (bottom) cross- 
correlations. 



and the finite length of the waveguide limited the number of bars that we could insert: 
we formed a scattering array of 40 bars in total. 

In contrast to [23], in the present case the bars were not only separated by 
random spacings d n (structural disorder), but also had random heights h n (compositional 
disorder). The bar heights spanned the range from 0.1 mm to 0.8 mm in steps of 0.1 mm. 
We derived sequences of bar heights h n and spacings d n by scaling the corresponding 
variables u n and A n of the Kronig-Penney model (JTj) which we used to obtain the 
transmission coefficient represented in figure |2] (see section [2]). In this way we obtained 
two random successions {h n } and {d n } with self-correlations of the form (JH]) with 
Ki = 0.207T and n 2 = 0.217T. This choice implied a rather high value of the transmittivity 
for most frequencies, with a sharp drop in the two narrow gaps corresponding to the 
intervals [0.20vr, 0.21tt] and [0.79tt, 0.80tt] of the left half of the first Brillouin zone. As 
in the example discussed in section [2], in addition to long-range self-correlations, the 
random sequences of bar heights and spacings were cross-correlated: we considered the 
cases of maximally positive (rj = vr/4) and negative (rj = — 7r/4) cross-correlations as 
well as the case in which heights and spacings of the bars were not cross-correlated 
( V = 0). 

In figure H] we show the arrangements of the scattering bars within the waveguide 
for the three cases discussed. The influence of the cross-correlations is clearly visible. In 
fact, positive cross-correlations tend to produce large spacings between high bars, while 
negative cross-correlation reduce the distance separating high bars. 

3.3. Transmission bands and Block vector calculation 

The measured transmission for the first mode versus the longitudinal wavevector is 
shown in figure [5j One can clearly see the first band. The band edges are the same 
for the cases of positive and absent cross-correlations, but in the case of negative 
cross-correlations the band is broader and the second gap is shifted towards higher 
wavenumbers k x . Note that in all the measurements we used the same bars, which 
were only shuffled and put in different positions in the three cases. This means that, 
although the constituents of the scattering arrangements are the same, cross-correlations 
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k x [1/m] 

Figure 5. (Colour online) Transmission coefficient T\ for the first propagating mode 
versus wave vector k x . The dotted line (blue) and the dashed line (red) respectively 
correspond to the cases of maximal positive and negative cross-correlations. The solid 
line (black) corresponds to the case in which no cross-correlations are present. 

can broaden the transmission band. We attribute this effect to a shadowing effect. In 
the case of negative cross-correlations, in fact, high bars are separated by much smaller 
distances than in the other cases. Therefore the waves do not "feel" the same underlying 
lattice structure in the case of negative cross- correlations and this is probably the origin 
of the enlarged band. 

To analyse the changes in the transmission gaps brought about by cross- correlations, 
we need to know the band structure of the experimental setup or, equivalently, the 
relation between the longitudinal wavevector k x and the Bloch vector k. 

We followed an empirical approach to establish a correspondence between k x and 
k. Specifically, we associated the edges of the transmission band to the values k = 
and k = it of the Bloch vector. We then used a linear interpolation to extend the k x — n 
correspondence to the whole Brilllouin zone. This corresponds to a linear fit k — ak, 
which is consistent with the fact that, as shown in [23], the finite width of the bars can 
be theoretically described by introducing an effective refractive index a so that k — > ak. 

4. Experimental results 

The theoretical results for the Kronig-Penney model, represented by the inverse 
localisation length ( TTTT) and the transmission coefficient pictured in figure [2j suggest 
that the transmission gap in the right neighbourhood of k = ir/5 should be reduced or 
enhanced according to whether the height-position correlations are positive or negative, 
whereas the contrary should happen for the gap in the left neighbourhood of k = Air/5. 
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Figure 6. (Colour online) Transmission coefficient T\ for the first propagating 
mode versus Bloch wavevector k. The dotted line (blue) and the dashed line 
(red) respectively correspond to the cases of maximal positive and negative cross- 
correlations. The shaded area highlights the variation of T\ due to cross-correlations. 

These expectations are confirmed by the experimental results, represented in figure El 
figure d and figure El 

In figure El we show the transmission coefficient T\ for the first propagating mode 
in the two extreme cases of maximum positive and negative cross-correlations between 
compositional and structural disorders. Figure [7] and figure |8] respectively compare 
the transmission coefficient T\ in the absence of cross-correlations (n = 0) with the 
same quantity measured when cross-correlations take the maximal positive (rj = 7r/4) 
and negative (77 = — it/ 4) values. The experimental data show that the transmission 
coefficient in the opaque gaps is reduced or enhanced according to the theoretical 
expectations. This is made particularly evident by the comparison of figure |2] and 
figure El In both cases switching from positive to negative cross-correlations decreases 
the transmittivity in the first opaque window and enhances it in the second gap. 
Figure [7] and figure [HI on the other hand, agree well with the behaviour of the inverse 
localisation length in the Kronig- Penney model represented in figure [U one can easily 
see, for instance, that positive cross-correlations decrease localisation in the first region 
of localised states and enhance it in the second interval while, correspondingly, they 
make the first transmission gap more shallow and deepen the second one. Negative 
cross-correlations have the same effects, but with reversed gaps. 
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Figure 7. (Colour online) Transmission coefficient T\ for the first propagating mode 
versus Bloch wavevector k. The clotted line (blue) and the solid line (black) respectively 
correspond to the cases of positive and absent cross-correlations. The shaded area 
highlights the variation of T\ due to positive cross-correlations. 




Figure 8. (Colour online) Transmission coefficient T\ for the first propagating mode 
versus Bloch wavevector k. The dashed line (red) and the solid line (black) respectively 
correspond to the cases of negative and absent cross-correlations. The shaded area 
highlights the variation of T\ due to negative cross-correlations. 
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5. Conclusions 

This work is the first experimental study of the effects that cross-correlations between 
different random potentials have on the transport properties of a ID waveguide. We 
analysed the microwave transmission of the first mode in a multimode waveguide with 
compositional and structural disorder where no mode coupling was present. Specifically, 
we inserted brass bars in the waveguide which acted as scatterers and were characterised 
by self- and cross-correlated random heights and positions. We showed how cross- 
correlations between the random heights and spacings of the scatterers can modulate 
the transmittivity in the non-transparent frequency windows created by long-range self- 
correlations of the disorder. The results agree with the theoretical predictions previously 
obtained for an aperiodic Kronig-Penney model. 

The agreement between theoretical predictions and experimental data is far 
from trivial, because the theoretical analysis was based on a number of simplifying 
assumptions, such as weak disorder and irrelevance of some experimental imperfections, 
e.g. absorption, which are present in any actual waveguide. Besides, the theoretical 
results for the transmission coefficient did not guarantee that the effect of cross- 
correlations would be visible also in ID devices of reduced length. The results obtained 
in the present experiments, therefore, support the validity of the theoretical predictions 
and show that the effects of cross-correlations are not eliminated by the reduced length 
of the random waveguide. 
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